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Experimental  investigation  of  bio-ethanol  autothermal  reforming  (ATR)  and  water-gas  shift  (WGS)  pro¬ 
cesses  for  hydrogen  production  and  regression  analysis  of  the  data  is  performed  in  the  study.  The  main 
goal  was  to  obtain  regression  relations  between  the  most  critical  dependent  variables  such  as  hydrogen, 
carbon  monoxide  and  methane  content  in  the  reformate  gas  and  independent  factors  such  as  air-to-fuel 
ratio  (A),  steam-to-carbon  ratio  (S/C),  inlet  temperature  of  reactants  into  reforming  process  (TATRin),  pres¬ 
sure  (p)  and  temperature  (Tatr)  in  the  ATR  reactor  from  the  experimental  data.  Purpose  of  the  regression 
models  is  to  provide  optimum  values  of  the  process  factors  that  give  the  maximum  amount  of  hydrogen. 
The  experimental  ATR  system  consisted  of  an  evaporator,  an  ATR  reactor  and  a  one-stage  WGS  reactor. 
Empirical  relations  between  hydrogen,  carbon  monoxide,  methane  content  and  the  controlling  parame¬ 
ters  downstream  of  the  ATR  reactor  are  shown  in  the  work.  The  optimization  results  show  that  within  the 
considered  range  of  the  process  factors  the  maximum  hydrogen  concentration  of  42  dry  vol.  %  and  yield 
of  3.8  mol  mol-1  of  ethanol  downstream  of  the  ATR  reactor  can  be  achieved  at  S/C  =  2.5,  A.  =  0.20-0.23, 
p  =  0.4  bar,  TATRin  =  230  °C,  Tatr  =  640  °C. 

©  2009  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  fuel  cell  technology  nowadays  is  one  of  the  most  promising 
high  efficiency  non-pollution  or  low-pollution  energy  supply  tech¬ 
nology,  and  can  play  an  important  role  in  promotion  of  distributed 
generation  of  electricity  as  well  as  cogeneration  of  electricity  and 
heat  [1,2].  Since  hydrogen  does  not  exist  on  Earth  in  pure  form, 
efficient  methods  of  producing  hydrogen  for  use  in  fuel  cell  tech¬ 
nologies  are  very  important  research  topics  [2,3].  Reforming  of 
hydrocarbons  and  alcohols  for  hydrogen  production  is  widely  dis¬ 
cussed  [4,3]  in  the  last  decade  since  these  sources  are  readily 
available  and  can  be  supplied  via  existing  infrastructure  [5].  Bio¬ 
ethanol  can  be  considered  as  a  renewable  and  basically  non-toxic 
product  that  may  be  produced  from  feedstock  and  therefore  is  one 
of  the  most  promising  sources  for  hydrogen  production  [2]. 

Generation  of  hydrogen-rich  gas  for  fuel  cells  from  alcohols  can 
be  done  using  three  main  reforming  processes:  steam  reforming 
(SR),  partial  oxidation  (POX)  and  autothermal  reforming  (ATR)  [4]. 
An  external  heat  source  is  required  in  the  SR  process  as  the  pro¬ 
cess  is  endothermic  but  it  gives  a  higher  hydrogen  concentration 
in  reformate  gas  compared  to  the  ATR  process  [1,6].  ATR  does  not 
require  an  external  heat  source  since  it  combines  endothermic  SR 
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and  exothermic  POX  processes  and  therefore  results  in  simpler 
design  and  higher  reforming  efficiencies  [1,6]. 

The  main  research  efforts  in  the  field  of  ATR  are  directed  at  devel¬ 
oping  and  optimizing  ethanol  ATR  catalysts,  the  fuel  conversion 
rates,  hydrogen  selectivity  and  materials  used  for  manufacturing 
[7-9].  There  are  a  number  of  studies  on  identifying  optimal  ATR 
reforming  conditions  with  no  determined  catalyst.  Ahmed  and 
Krumpelt  [6]  theoretically  study  the  reforming  reaction  and  have 
found  that  the  maximum  theoretical  reforming  efficiency  can  be 
achieved  at  the  thermoneutral  point,  at  which  for  case  of  ethanol 
ATR  the  value  of  A  is  0.2  and  S/C  =  0.9.  Semelsberger  et  al.  [10]  have 
also  theoretically  estimated  the  optimal  conditions  for  maximiz¬ 
ing  hydrogen  concentration  downstream  of  the  ATR  reactor  with 
ethanol  as  a  fuel.  They  have  done  simulation  at  thermoneutral  and 
equilibrium  conditions  and  the  results  showed  that  to  achieve  a 
hydrogen  concentration  of  41-43%  in  the  reformate  gas  the  reac¬ 
tion  temperature  should  be  in  the  range  of  257-327  °C  and  S/C  from 
1.6  to  2.9.  Rabenstein  and  Hacker  [11]  have  analyzed  the  hydro¬ 
gen  yields  and  total  energy  demand  for  the  process,  as  well  as 
determined  the  coke  formation  region,  at  a  wide  range  of  S/C,  A 
and  reforming  temperatures  using  Gibb’s  free  energy  minimization 
thermodynamical  calculations.  Hagh  [12]  has  estimated  that  for 
nearly  maximum  thermal  efficiency  of  the  system  it  should  operate 
at  inlet  temperature  of  reactants  in  the  ATR  reactor  TATRin  around 
670  °C,  A  =  0.325  and  S/C  =  3.  Ersoz  et  al.  [  13 ]  simulated  an  adiabatic 
ethanol  autothermal  reforming  process  at  chemical  equilibrium 
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Nomenclature 

ATR 

autothermal  reforming 

ch4ATR 

methane  concentration  in  reformate  gas  down¬ 
stream  of  ATR  (%  vol.dry  basis) 

CH4WGs 

methane  concentration  in  reformate  gas  down¬ 
stream  of  WGS  (%  vol.dry  basis) 

COatr 

carbon  monoxide  concentration  in  reformate  gas 
downstream  of  ATR  (%  vol.dry  basis) 

COwgs 

carbon  monoxide  concentration  in  reformate  gas 
downstream  of  WGS  (%  vol.dry  basis) 

ghsvatr 

gas  hourly  space  velocity  in  ATR  reactor  (h-1 ) 

GHSVwgs  gas  hourly  space  velocity  in  WGS  reactor  (h-1 ) 

H2ATR 

hydrogen  concentration  (%  vol.dry  basis)  or  yield 
(mol  of  hydrogen  produced  per  mol  of  ethanol  sup¬ 
plied  to  reactor)  in  reformate  gas  downstream  of 
ATR 

H2WGS 

hydrogen  concentration  in  reformate  gas  down¬ 
stream  of  WGS  (%  vol.dry  basis) 

A  Hr 

heat  of  reaction  (kj  mol-1 ) 

V 

system  gauge  pressure  (bar) 

PEMFC 

proton  exchange  membrane  fuel  cell 

S/C 

steam-to-carbon  feed  ratio  (mol  mol-1 ) 

SR 

steam  reforming 

Tatr 

average  temperature  of  reactants  in  the  catalyst  of 
ATR  reactor  (°C) 

^ATRin 

inlet  temperature  of  reactants  into  ATR  reactor  (°C) 

^ATRout 

temperature  of  reactants  at  the  exit  of  the  catalyst 
of  ATR  reactor  (°C) 

^WGSin 

temperature  of  reactants  at  the  entrance  of  the  cat¬ 
alyst  of  WGS  reactor  (°C) 

fwGSout 

temperature  of  reactants  at  the  exit  of  the  catalyst 
of  WGS  reactor  (°C) 

WGS 

water-gas  shift 

VIF 

variance  inflation  factor 

Greek  letters 

X 

air-to-fuel  combustion  ratio 

conditions,  and  have  estimated  that  the  optimal  conditions  of  the 
reforming  process  for  maximizing  the  hydrogen  concentration  and 
minimizing  carbon  monoxide  concentration  and  additional  heat 
supply  to  the  system  can  be  achieved  at  S/C  value  of  3.5,  Tatr  value 
of  700  °C  and  oxygen-to-carbon  ratio  of  0.48.  The  achieved  maxi¬ 
mum  values  of  hydrogen  and  thermal  efficiency  of  the  system  were 
32%  and  78%  respectively. 

As  can  be  seen,  the  results  considerably  differ  since  the  con¬ 
sidered  constraints  and  optimal  parameters  in  each  research  work 
were  different.  Empirical  models  derived  from  experimental  results 
can  be  used  for  detailed  studies  and  process  simulation.  The 
advantage  of  such  models  is  that  they  can  predict  and  show  the 
interactions  between  all  variables  of  the  reforming  system  which 
is  difficult  or  even  impossible  to  describe  with  mechanistic  models 
[14].  There  are  many  studies  performed  in  the  field  of  modeling  and 
optimization  of  chemical  reforming  processes  by  using  empirical 
models  [15-17].  Larentis  et  al.  [16]  have  compared  two  mathemat¬ 
ical  modeling  approaches:  empirical  and  phenomenological.  The 
empirical  model  was  found  to  be  more  efficient,  simpler  and  led 
to  better  results  than  those  obtained  with  the  phenomenological 
model  approach  in  modeling  and  optimization  of  reforming  pro¬ 
cesses.  However,  currently  there  is  a  lack  of  experimental  data  and 
empirical  models  for  modeling  and  optimization  of  the  bio-ethanol 
ATR  reforming  process. 


The  objective  of  the  present  work  was  to  obtain  regres¬ 
sion  relations  between  the  response  variables  of  interest  which 
are— hydrogen,  carbon  monoxide  and  methane  content  in  the  refor¬ 
mate  gas  downstream  of  the  ATR  reactor  and  the  most  critical 
process  factors  such  as  S/C,  X,  p,  TATRin,  Tatr  and  GHSVATR  by  using 
data  from  experimental  investigations  of  the  bio-ethanol  reform¬ 
ing  process.  The  purpose  of  the  regression  models  was  to  determine 
the  optimal  values  of  the  process  factors  which  maximize  hydrogen 
concentration  and  yield  (defined  as  the  molar  ratio  of  the  obtained 
hydrogen  to  the  ethanol  feed)  downstream  of  the  ATR  reactor.  The 
experiments  were  made  at  Fraunhofer  Institute  for  Solar  Energy 
systems  (ISE). 

2.  Experimental 

2.1.  Catalyst 

Commercial  honeycomb  catalysts  were  used  in  the  ATR  reactor 
(with  the  optimum  temperature  of  operation  around  700  °C)  and 
in  the  medium-temperature  WGS  reactor  (with  the  optimal  inlet 
temperature  of  reactants  in  the  range  of  250-260  °C). 

2.2.  Apparatus 

The  autothermal  reforming  of  bio-ethanol  over  the  catalyst  was 
carried  out  in  a  cylindrical  reactor  with  cross-section  diameter  of 
26.0  mm  built  of  steel  with  high  heat  resistance.  The  catalyst  bed 
with  diameter  of  20.6  mm  and  height  of  41.2  mm  was  placed  in 
the  reactor  as  shown  in  Fig.  1.  The  WGS  catalyst  with  dimensions 
77.1  mm  x  59.0  mm  x  18.3  mm  was  placed  in  the  middle  of  the  WGS 
reactor  as  shown  in  Fig.  2. 

When  starting  the  experiments,  the  ATR  and  WGS  reactor  cat¬ 
alysts  were  pre-heated  in  a  flow  of  nitrogen  until  the  catalysts 
reached  the  temperature  above  100  °C.  The  ATR  catalyst  was  then 


Fig.  1.  Schematic  view  of  the  ATR  reactor  and  positioning  of  the  catalyst  (dimensions 
are  given  in  mm)  [18]. 
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Fig.  3.  Results  of  measurements  of  the  component  concentrations  downstream  of  ATR  and  WGS  reactors  at  S/C  =  1.5,  A.  =  0.3,  Tatrouc  ^700°C,  rWGSin^345°C  and 
Twcsout  S»400°C.  Solid  lines  present  the  values  measured  with  gas  analyzer  and  symbols— the  values  measured  with  gas  chromatograph. 


heated  with  water  vapor  flow  to  circa  220  °C,  and  the  reactants  were 
subsequently  introduced  into  the  reactor.  After  the  ATR  reactor  had 
reached  stable  temperature  and  pressure  in  the  system  the  refor¬ 
mate  gas  (reaction  products  of  the  ATR  reactor)  was  directed  into 
the  WGS  reactor.  As  the  water-gas  shift  reaction  is  only  slightly 
exothermic  the  additional  reactor  heating  was  used  in  order  to 
compensate  for  the  heat  losses  to  surroundings  and  maintain  the 
required  temperature  in  the  reactor. 

The  reactant  feed  was  a  mixture  of  distilled  water  and  bio¬ 
ethanol  (99%  of  bio-ethanol;  1%  of  methylethylketone;  10  ppm  of 
denatonium  benzoate)  which  was  pre-heated  and  vaporized  in  the 
vaporizer,  and  air,  pre-heated  in  the  tube  bent  around  the  vapor¬ 
izer.  Temperatures  of  the  reactants  and  the  catalyst  were  measured 
with  “type-K”  thermocouples,  which  were  placed  at  the  inlet  of 
ATR  reactor  and  in  the  catalyst  in  three  locations:  at  the  inlet,  in 
the  middle  and  at  the  exit  of  the  catalyst.  The  average  of  these 
temperatures  was  calculated  and  used  as  the  factor  representing 
temperature  of  the  ATR  catalyst  for  deriving  of  multiple  regression 
models  from  the  experimental  data.  Thermocouples  were  placed 


also  in  two  locations  on  the  WGS  reactor  catalyst:  at  the  inlet  and 
at  the  exit  of  the  catalyst  (Fig.  2).  The  airflow  was  controlled  with 
a  mass  flow  meter,  model  “1259X-20000V”  (“MKS  Instruments”). 
The  pump  “HPLH20/200  PF”  (“CAT”)  was  used  for  water  and  ethanol 
mixture  dosing.  The  product  gas  composition  was  measured  by  an 
online  gas  analyzer,  model  “Optima  A02000”  (“ABB”)  and  a  gas 
chromatograph,  model  “Agilent  6890  Series  LC/MSD”  (“Agilent”). 
The  following  components  were  measured  in  the  reaction  products 
downstream  of  the  ATR  and  WGS  reactors:  hydrogen  (H2),  car¬ 
bon  monoxide  (CO),  carbon  dioxide  (C02),  methane  (CH4),  nitrogen 
(N2),  oxygen  (02)  and  water. 

2.3.  Experimental  conditions 

Experimental  conditions  were  selected  based  on  the  previous 
experience  in  research  of  ATR  processes  at  Fraunhofer  Institute  for 
Solar  Energy  systems  and  results  of  chemical  equilibrium  simula¬ 
tions  done  with  software  “ChemCAD”,  where  the  optimal  S/C  and 
A  ranges  according  to  the  ATR  catalyst  optimum  temperature  of 
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Table  1 

The  process  factors  and  response  variables. 


The  process  factors 

Range 

o 

> 

S/C 

0.5-2. 5 

■a 

k 

0.19-0.36 

c 

o 

P 

0.12-0.46 

55 

TvTRin 

230-360 

o 

Q. 

Taj  r 

641-782 

E 

GHSVAtr 

14830-33341 

o 

o 

T/VGSin 

360-380 

w 

TO 

fwGSout 

265-435 

O 

GHSVwgs 

2405-4164 

The  response  variables:  H2atr,  COAtr,  CH4Atr. 


operation  were  estimated.  The  present  ranges  of  process  factors 
were  chosen  sufficiently  broad  to  obtain  optimal  process  con¬ 
ditions.  The  range  of  process  factors  and  the  chosen  response 
variables  are  given  in  Table  1 .  The  developed  ATR  and  WGS  reac¬ 
tor  system  was  a  part  of  a  research  project  in  which  a  bio-ethanol 
fuel  processing  system  was  designed  to  be  coupled  with  a  proton 
exchange  membrane  fuel  cell  (PEMFC).  The  aim  of  these  exper¬ 
iments  was  to  achieve  a  stable  reforming  process  in  which  the 
measured  gas  concentrations  could  be  compared  to  the  results 
obtained  from  chemical  equilibrium  calculations,  and  to  minimize 
the  amount  of  additional  heating  of  the  system,  and  achieve  max¬ 
imum  hydrogen  and  minimum  carbon  monoxide  concentrations 
in  the  reformat  gas  downstream  of  the  ATR  and  WGS  reactors  as 
PEMFC  is  highly  sensitive  to  carbon  monoxide  [1  ]. 

3.  Results  and  discussions 

3.1  Experimental  results 

Experimental  results  obtained  during  one  of  the  experimental 
runs  are  shown  in  Fig.  3.  The  experiment  consisted  of  coupled  ATR 
and  WGS  reactors.  After  coupling  with  the  WGS  reactor  quite  rapid 
stabilization  of  gas  concentrations  was  achieved.  Measurements  of 
the  concentrations  of  components  (Fig.  3)  downstream  of  the  ATR 
reactor  showed  that  37%  of  H2,  12%  CO  and  approximately  1%  of 
CH4  concentrations  can  be  obtained  which  agree  quite  well  with 
theoretical  calculations  [  10,13  ].  Connection  of  a  WGS  reactor  down¬ 
stream  of  the  ATR  allows  reducing  CO  concentration  to  about  3.5% 
(Fig.  3)  and  increasing  H2  concentration  to  40%  as  well  as  raising 
C02  concentration  from  13%  to  19%  according  to  the  water-gas  shift 


Fig.  4.  Comparison  of  gas  concentrations  obtained  from  chemical  equilibrium  (eq) 
calculations  and  experiments  (exp)  (at  k  =  0.29;  TATRout  =  700  °C;  p  =  1.26  bar). 

reaction  (Eq.  (1)). 

C0  +  H20^C02  +  H2,  AHr  =  -41  kjmor1  (1) 

Additionally,  CH4  concentration  increases  to  about  2%  down¬ 
stream  of  WGS  reactor  indicating  that  other  reactions  might  be 
occurring  at  the  same  time,  and  one  of  the  possible  reactions  is 
the  following: 

2CO  +  2H2  ^CH4  +  C02,  AHr  =  -247.1  kjmor1  (2) 

The  obtained  CO  concentration  downstream  of  the  WGS  is  in 
quite  a  good  agreement  with  other  WGS  reactor  studies  [19],  but 
does  not  satisfy  the  requirements  of  PEMFC  [1]  and  additional  CO 
removal  reactors  are  necessary.  Measurements  (Fig.  3)  show  that 
both  ATR  and  WGS  processes  are  stable  and  addition  of  the  WGS 
reactor  rapidly  changes  composition  of  the  reformate  gas  which  is 
important  for  dynamic  change  of  fuel  cell  capacity. 

Comparison  of  the  experimental  results  with  results  of  chemical 
equilibrium  calculations  for  different  S/C  ratios  (Fig.  4)  show  that 
the  experimental  results  are  in  fairly  good  agreement  with  chemical 
equilibrium  calculations.  Gibb’s  free  enthalpy  equilibrium  reactor 
and  isothermal  conditions  were  chosen  for  the  simulation. 

3.2.  Statistical  modeling 

Before  the  regression  analysis  rough  errors  in  the  measured 
results  were  removed  and  it  was  checked  that  the  obtained  results 
of  the  response  variables  correspond  to  normal  distribution  since 


Table  2 

Correlation  coefficients  between  hydrogen  concentration  (%  vol.dry  basis)  in  the  reaction  products  measured  downstream  of  the  ATR  reactor  and  process  factors  as  well 
as  between  process  factors  (first  number  is  the  value  of  correlation  coefficient;  the  number  in  parenthesis  is  p-value,  indicating  statistical  significance  of  the  estimated 
correlations). 


k  p 

GHSVatr 

S/C 

TvTRin 

T/ \TR 

H2ATR 

-0.45  (<0.01 )  0.67  (<0.01 ) 

0.54  (<0.01 ) 

-0.21  (0.01) 

-0.65  (<0.01 ) 

-0.33  (<0.01 ) 

k 

-0.30  (<0.01 ) 

-0.63  (<0.01 ) 

0.77  (<0.01 ) 

0.28  (<0.01 ) 

0.08  (0.36) 

P 

GHSVAtr 

S/C 

TvTRin 

0.46  (<0.01 ) 

-0.24(0.003) 

-0.74  (<0.01 ) 

-0.44  (<0.01 ) 

-0.32  (<0.01 ) 

<0.01  (0.98) 

-0.28  (<0.01 ) 
0.21  (<0.01 ) 
-0.30  (<0.01 ) 
0.52  (<0.01 ) 

Table  3 

Correlation  coefficients  between  carbon  monoxide  concentration  (%  vol.dry  basis)  in  the  reaction  products  measured  downstream  of  the  ATR  reactor  and  process  factors  and 
correlation  coefficients  between  the  process  factors. 


A.  p 

GHSVatr 

S/C 

^ATRin 

Tatr 

COATR  -0.62  (<0.01 )  -0.06  (0.50) 

k  -0.35  (<0.01 ) 

P 

GHSVatr 

S/C 

TvTRin 

0.64  (<0.01 ) 

-0.69  (<0.01 ) 

0.49  (<0.01 ) 

-0.86  (<0.01 ) 

0.79  (<0.01 ) 

-0.25  (<0.01 ) 

-0.74  (<0.01 ) 

0.27  (<0.01 ) 

0.240  (<0.01 ) 

-0.46  (<0.01 ) 

-0.42  (<0.01 ) 

0.02  (0.77) 

0.54  (<0.01 ) 
-0.09  (0.26) 
-0.33  (<0.01 ) 
0.21  (<0.01 ) 
-0.44  (<0.01 ) 
0.31  (<0.01 ) 
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Table  4 

Correlation  coefficients  between  methane  concentration  (%  vol.dry  basis)  in  the  reaction  products  measured  downstream  of  the  ATR  reactor  and  process  factors  and  correlation 
coefficients  between  the  process  factors. 


A.  p 

GHSVatr 

S/C 

TvTRin 

Tatr 

CH4Atr 

A 

P 

GHSVatr 

S/C 

TvTRin 

-0.30  (<0.01 )  -0.26  (<0.01 ) 

-0.35  (<0.01 ) 

0.18  (0.04) 

-0.63  (<0.01 ) 

0.52  (<0.01 ) 

-0.50  (<0.01 ) 

0.78  (<0.01 ) 

-0.29  (<0.01 ) 

-0.74  (<0.01 ) 

0.51  (<0.01 ) 

0.34  (<0.01 ) 

-0.44  (<0.01 ) 

-0.37  (<0.01 ) 

0.05  (0.59) 

0.28  (<0.01 ) 
0.17  (0.05) 
-0.26  (<0.01 ) 
0.20(0.02) 
-0.27  (<0.01 ) 
0.44  (<0.01 ) 

Table  5 

Parameter  estimates  of  the  regression  of  hydrogen  concentration  (H2atr)  (the  model 
is  fitted  with  146  data  points;  p-values  <0.01  for  all  terms  therefore  are  not  shown 
in  the  table). 


Term 

Estimate 

Standard  error 

t- ratio 

Constant 

-72.43 

8.18 

-8.85 

GHSVatr  x  TATRin 

-39.25  x  10“7 

5.47  x  10“7 

-7.18 

GHSVatr  x  TAtr 

19.86  x  10“7 

2.34  x  10“7 

8.47 

A  x  A 

-124.33 

13.37 

-9.30 

Ax  S/C 

23.85 

3.30 

7.24 

P 

22.20 

3.84 

5.77 

P  x  GHSVatr 

-85.30  xl0“5 

14.56  x  10“5 

-5.86 

S/C  x  S/C 

-1.86 

0.26 

-7.12 

TvTRin 

0.85 

0.06 

13.27 

TvTRin  x  rATRin 

-93.13  x  10-5 

7.43  x  10“5 

-12.53 

7/VTRin  X  TAtr 

-28.71  x  10“5 

2.78  x  10-5 

-10.32 

that  is  the  condition  for  use  of  these  data  in  regression  analysis.  The 
statistics  program  “STATGRAPHICS  Plus  for  Windows  2.1  ”  was  used 
for  statistical  analysis  and  empirical  model  building.  Correlation 
coefficients  between  the  response  variables  and  the  process  factors 
used  in  the  regression  analysis  as  well  as  between  the  process  fac¬ 
tors  themselves  measuring  the  strength  of  the  linear  relationship 
between  these  variables,  obtained  for  the  results  of  ATR  reactor, 
are  shown  in  Tables  2-4.  The  linear  correlation  of  the  process  fac¬ 
tors  and  response  variables  are  statistically  significant  at  the  95% 
confidence  level  except  for  the  correlation  between  p  and  COATR. 
Chemical  equilibrium  calculations  show  that  within  the  studied 
range  of  ATR  reactor  temperatures  CO  content  in  the  reformate  gas 
is  increasing  relatively  fast,  presumably  mainly  due  to  the  water-gas 
shift  reaction.  However,  the  impact  of  pressure  on  the  equilibrium 
product  composition  of  water-gas  shift  reaction  within  this  temper¬ 
ature  range  is  very  small.  The  above  considerations  could  explain 
the  absence  of  linear  correlation  between  p  and  COATR  (Table  3). 
The  results  of  correlation  coefficients  also  show  that  collinearity 


Table  5a 

ANOVA  and  other  diagnostic  statistics  for  H2atr- 


Table  6 

The  parameter  estimates  of  the  regression  of  carbon  monoxide  concentration  (COatr)  (the  model  is  fitted  with  144  data  points;  p-values  <0.01  for  all  terms  therefore  are  not 
shown  in  the  table). 


Table  6a 

ANOVA  and  other  diagnostic  statistics  for  COatr. 
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Table  7 

The  parameter  estimates  of  the  regression  of  methane  concentration  (CH4Atr)  (the 
model  is  fitted  with  135  data  points;  p-values  <0.01  for  all  terms  therefore  are  not 
shown  in  the  table). 


Term 

Estimate 

Standard  error 

t-  ratio 

Constant 

29.33 

2.88 

10.18 

GHSVAtr  x  GHSVatr 

-1.74  x  10“9 

3  x  10-10 

-5.80 

X 

69.29 

12.50 

5.55 

X  X  GHSVatr 

37.15  x  10“5 

5.21  x  10“5 

7.21 

x  x  rATR 

-0.102 

0.017 

-6.004 

p 

-34.70 

5.67 

-6.12 

p  x  GHSVatr 

10.76  x  10“5 

2.29  x  10-5 

4.70 

V  x  Latr 

0.041 

0.007 

5.47 

S/C 

-3.56 

0.76 

-4.66 

S/C  X  S/C 

0.169 

0.026 

6.43 

S/C  x  p 

1.40 

0.231 

6.07 

S/C  x  GHSVatr 

-2.82  x  10“5 

40.03  x  10“7 

-7.03 

s/c  x  rATR 

0.004 

94.01  x  10“5 

4.63 

FvTRin 

-0.145 

0.014 

-10.05 

LATRin  x  f ATR  in 

7.81  x  10“5 

72.72  x  10“7 

10.74 

FvTRin  X  TATR 

13.59  x  10-5 

2.01  x  10“5 

6.77 

Latr 

-0.029 

0.004 

-6.96 

between  X  and  S/C,  X  and  GHSVATR,  S/C  and  GHSVATR  is  quite  strong 
which  means  that  the  interaction  terms  of  the  process  factors  are 
important  for  regression  analysis.  The  collinearity  between  the  pro¬ 
cess  factors  also  indicated  that  it  may  complicate  use  of  the  obtained 
regression  equations  in  optimization  calculations  without  correct¬ 
ing  measures,  since  individual  effects  of  multicollinear  factors  on 
the  response  variables  may  not  be  determined. 

It  was  found  from  chemical  equilibrium  calculations  that  the 
relations  between  most  of  the  factors  and  response  variables  are 
nonlinear  and  can  be  approximated  quite  well  with  the  second- 
order  polynomial  function  within  the  considered  range  of  values. 
During  the  process  of  selection  of  the  appropriate  regression  model, 
multiple  regression  models  containing  the  process  factors  in  lin¬ 
ear  form  only  were  compared  with  the  models  obtained  from 
full  second-order  polynomial  equations.  It  was  found  that  for  all 
response  variables  the  models  obtained  from  the  second-order 
polynomial  equations  provide  better  fit  in  terms  of  some  of  the 
main  statistical  parameters,  namely,  these  equations  had  higher 
values  of  coefficient  of  determination  adjusted  for  degrees  of  free¬ 
dom  (“ R 2  adjusted”)  and  F-ratios,  and  lower  values  of  mean  square 
errors  (MSE)  than  the  linear  regression  equations.  Since  authors  did 
not  have  any  theoretical  reason  for  using  another  nonlinear  form  of 
the  regression  model  equation  to  describe  response  surfaces,  the 
second-order  polynomial  equation  was  used  for  fitting  the  experi¬ 
mental  data  for  all  response  variables.  Low-order  polynomial  form 
of  relationship  is  usually  used  in  response  surface  methodology 
problems  [20]  and  successfully  employed  in  regression  analysis  of 
chemical  processes  [15-17,21].  The  procedure  started  with  using 
a  complete  form  of  the  second-order  model  including  all  linear, 
second-order  terms  and  all  two-factor  interaction  terms  and  per¬ 
forming  stepwise  backward  elimination  afterwards.  Statistically 
non-significant  terms  were  eliminated  from  the  model  according 
to  p-values  and  t- statistics  by  eliminating  the  terms  with  highest 
p-values  and  lowest  absolute  t- ratios  first.  Each  time  after  elimina¬ 
tion  of  the  term,  the  regression  equation  is  recalculated,  and  the 
elimination  procedure  continues  until  all  terms  left  in  the  regres¬ 
sion  model  are  statistically  significant  at  95%  confidence  level  (the 
significance  level  p  <  0.05). 


Table  8 

The  parameter  estimates  of  the  ridge  regression  of  hydrogen  concentration  (H2atr, 
%  vol.dry  basis)  (ridge  parameter  =  0.016,  MSE  =  0.74,  R 2  adjusted  for  degree  of  free¬ 
dom  =  69.7%). 


Results  of  the  reduced  regression  models  with  all  the  insignifi¬ 
cant  terms  eliminated  from  the  regression  model  for  the  response 
variables  H2AtR,  COatr  and  CH4ATR  are  presented  in  Tables  5-7, 
Tables  5a-7a. 

Values  of  determination  coefficients  “R2  adjusted”  show  that 
the  variation  in  response  variables  is  quite  well  explained  by  the 
regression  models.  Since  F-ratio  for  all  regression  models  exceed 
the  critical  values  at  the  chosen  significance  level,  it  can  be  con¬ 
cluded  that  the  models  adequately  describe  the  experimental  data. 
Since  the  Durbin-Watson  statistics  is  greater  than  1.4  for  all  models, 
it  can  be  concluded  that  there  is  probably  no  serious  autocorrela¬ 
tion  in  the  residuals.  Plots  of  the  residuals  versus  row  number  in 
which  they  appear  (Fig.  5)  indicated  that  there  was  no  clear  pattern 
and  that  distribution  of  the  residuals  was  approximately  normal. 
Each  of  the  observed  value  is  compared  with  the  predicted  value 
calculated  from  the  model,  as  shown  in  Fig.  6. 

3.3.  Ridge  regression 

Since  the  values  of  variance  inflation  factors  (VIF)  for  the  regres¬ 
sion  coefficients  of  the  obtained  regression  models  were  much 
higher  than  the  assumed  “cut-off’  value  of  10,  indicating  that  the 
multicollinearity  among  the  process  factors  may  cause  a  problem 
in  determination  of  individual  effects  of  the  process  factors  on  the 
response  variable,  ridge  regression  was  used  to  obtain  the  regres¬ 
sion  models  for  finding  the  optimum  values  of  process  factors  of 
ATR  process.  The  ridge  parameters  were  chosen  with  condition 
that  none  of  the  VIF  values  of  the  regression  coefficients  exceed  10. 
Tables  8  and  9  show  the  estimates  of  regression  coefficients  given 
by  the  ridge  regression  and  the  corresponding  VIF  values. 

3.4.  Optimization 

By  using  the  obtained  ridge  regressions  for  hydrogen  concen¬ 
tration  and  yield  downstream  of  the  ATR  reactor  the  optimal 
conditions  for  ATR  reforming  process  with  the  aim  of  getting 
maximum  hydrogen  content  in  the  reformat  gas  and  maximum 
hydrogen  yield  “MS  Excel  Solver”  tool  has  been  used.  The  range 
of  the  process  factors  which  are  covered  in  the  experiments  and 
define  the  validity  range  of  the  obtained  regression  models  were 
used  as  the  constraints  in  the  optimization  model  (Table  1 ). 


Table  7a 

ANOVA  and  other  diagnostic  statistics  for  CH4Atr. 
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Fig.  5.  Plots  of  residual  distributions  versus  row  numbers  for  regression  models  of  (a)  hydrogen  concentration,  (b)  carbon  monoxide  concentration,  and  (c)  methane 
concentration. 


Fig.  6.  Comparison  between  predicted  and  observed  values  of  (a)  hydrogen  concentration,  (b)  carbon  monoxide  concentration,  and  (c)  methane  concentration. 


Table  9 

The  parameter  estimates  of  the  ridge  regression  of  hydrogen  yield  (mol  H2atr 
mol-1  C2H5OH)  (ridge  parameter  =  0.015,  MSE  =  0.008,  R2  adjusted  for  degree  of  free¬ 
dom  =76.1%). 


The  obtained  results  of  the  maximum  hydrogen  concentration 
in  the  reforming  products  and  yield  as  well  as  values  of  the  cor¬ 
responding  optimum  process  parameters  and  thermal  efficiency 
(calculated  as  the  ratio  of  the  heat  of  hydrogen  obtained  to  the  heat 
of  ethanol  feed)  are  presented  in  Table  10.  The  maximum  hydrogen 
yield  at  lower  A  values  is  expected  from  thermodynamic  analysis 
therefore  it  is  not  surprising  to  see  that  the  optimum  value  of  A 
(Table  10)  corresponds  to  the  minimum  A  value  within  the  studied 


Table  10 

Summary  of  the  optimal  operational  conditions  for  bio-ethanol  ATR  process. 


Performance  criteria 

Maximum  H2atr  concentration 
(%  vol.dry  basis)  in  the  products 

Maximum  H2Atr  yield  (mol 
H2ATR  mol-'CjHsOH) 

Thermal 
efficiency  (%) 

42 

3.8 

74 

Optimal  operational  parameters 

S/C  A 

P  TATRin  (°C) 

Tatr  (°C) 

2.5  0.20-0.23 

0.4  230 

640 

range.  It  is  also  expected  theoretically  that  increase  of  S/C  will  result 
in  higher  hydrogen  yields.  Therefore,  the  optimum  value  of  S/C  is 
obtained  at  the  maximum  value  in  the  studied  range  of  parame¬ 
ters.  The  results  of  the  optimum  values  of  Tatr  and  S/C  confirm  the 
shift  of  maximum  hydrogen  yield  towards  lower  values  of  Tatr  at 
the  higher  S/C  values  which  is  also  seen  from  chemical  equilibrium 
calculations  performed  by  the  authors  and  other  researchers  [11]. 
Since  lower  Tatr  values  correspond  to  lower  temperatures  of  reac¬ 
tants  at  the  inlet  of  ATR  reactor  ^ATRin*  the  optimum  value  of  TATRin 
also  equals  to  the  lowest  in  the  studied  range.  The  effect  of  pressure 
on  the  hydrogen  yield  contradicts  to  the  expectation  since  higher 
values  of  hydrogen  yield  are  expected  at  lower  values  of  pressure 
according  to  the  thermodynamic  analysis.  One  of  the  reasons  for 
the  contradiction  could  be  the  possibility  that  individual  effect  of 
pressure  on  the  hydrogen  yield  is  not  correctly  clarified  due  to  the 
experimental  conditions. 

For  operation  of  real  ATR  systems  the  obtained  optimum  values 
of  A,  S/C  and  ATR  temperatures  mean  that  operation  conditions  at 
the  maximum  hydrogen  yield  correspond  to  the  range  where  coke 
formation  is  not  expected  according  to  thermodynamic  analysis. 
Higher  S/C  values  mean  that  more  heat  is  needed  for  heating  and 
evaporation  of  the  supplied  water  to  the  ATR  reactor.  Therefore, 
heat  balance  of  the  ATR  reactor  has  to  be  considered  in  order  to 
determine  the  optimum  operating  conditions  of  reactor  from  the 
energy  efficiency  viewpoint. 

As  the  maximum  theoretical  yield  of  hydrogen  which  can  be 
obtained  for  ATR  process  is  4.8  mol  mol-1  of  ethanol  feed  and  the 
maximum  theoretical  ethanol  ATR  reformer  efficiency  at  the  ther¬ 
moneutral  conditions  is  93.7%  [6]  there  is  a  possibility  for  further 
improvement  of  the  ATR  process.  Apart  from  the  search  of  more  effi¬ 
cient  catalysts  which  would  improve  composition  of  the  reformed 
gas,  optimization  of  the  ATR  and  subsequent  WGS  reactor  processes 
can  be  done  by  experimental  investigation  of  optimal  values  of  S/C, 
A,  temperatures  and  pressure  both  with  respect  to  composition  of 
the  reforming  products  and  the  overall  reforming  system  efficiency. 
A  well  designed  experimental  work  to  find  the  optimal  values  of 
S/C  within  the  range  of  1.5-3,  A  within  the  range  of  0.2-0.3  and 
temperature  of  ATR  within  range  of  600-700  °C  should  be  done  in 
the  future  to  optimize  performance  of  the  ATR  reactor  considering 
also  minimization  of  CO  and  CH4  content  in  the  product  gas  in  the 
objective  function. 
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4.  Conclusions 

The  experimental  results  of  bio-ethanol  autothermal  reform¬ 
ing  process  are  in  good  agreement  with  the  results  of  chemical 
equilibrium  calculations.  Multiple  regression  analysis  of  the  mea¬ 
sured  reformate  gas  composition  data  downstream  of  the  ATR 
reactor  shows  that  regression  models  with  a  reasonable  fit  to  the 
experimental  data  can  be  obtained.  These  models  can  be  used 
for  evaluation  of  H2,  CO  and  CH4  concentrations  in  the  reformate 
gas,  depending  on  the  values  of  S/C,  X,  p,  TATRin,  Tatr  and  GHSV 
within  the  range  of  values  of  regressors.  Optimization  performed 
with  the  obtained  regression  equations  shows  that  the  maximum 
H2ATR  concentration  downstream  of  ATR  reactor  of  42%  and  maxi¬ 
mum  H2AtR  yield  of  3.8  can  be  achieved  at  S/C  =  2.5,  X  =  0.20-0.23, 
^ATRin  =  230  °C,  TAxr  =  640  °C. 
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